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f^ Abstract We study a model for the translocation of proteins across mem- 

CN branes through a nanopore using a ratcheting mechanism. When the protein 

1— H enters the nanopore it diffuses in and out of the pore according to a Brownian 

, ^ motion. Moreover, it is bound by ratcheting molecules which hinder the diffu- 

sion of the protein out of the nanopore, i.e. the Brownian motion is reflected 
such that no ratcheting molecule exits the pore. New ratcheting molecules bind 



at rate 7. Extending our previous approach (Depperschmidt and Pfaffelhuber 



r— I 2010) we allow the ratcheting molecules to dissociate (at rate 6) from the 

Pm protein (Model I). We also provide an approximate model (Model II) which 

Ph assumes a Poisson equilibrium of ratcheting molecules on one side of the cur- 

(~| rent reflection boundary. Using analytical methods and simulations we show 

that the speed of both models are approximately the same. Our analytical 

results on Model II give the speed of translocation by means of a solution of 



an ordinary differential equation. 
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1 Introduction 



Most proteins are generated within the cellular cytosol and have to be trans- 
ported where they are needed. Frequently, they have to be transported across 
membranes, e.g. into the endoplasmatic reticulum, mitochondria or other or- 
ganelles, or they are secreted, i.e. transported across the cell wall (jNeupert 
and Brunner 2002 1. Protein translocation occurs either co-translationally, e.g. 



when ribosomes attach to the membrane of the endoplasmatic reticulum and 
proteins are transported into the lumen through a nanopore already during 



translation, or post-translationally (Rapoport 2007). 



Consider translocation into mitochondria first. Over 99% of mitochondrial 



proteins are translocated post-translationally (Wickner and Schekman 2005) 



and the molecular mechanisms have been studied in detail. When a substrate 
is translocated into the mitochondrium, it has to cross the translocase outer 
membrane (TOM) and the translocase inner membrane (TIM). The mito- 
chondrial heat shock protein 70 (mtHsp70) is known to play a crucial role in 
translocation of substrates ( Click 1995[ ) . mtHsp70 changes between its ATP- 
bound form, where it has an open pocket for binding to the substrate, and 



the ADP-bound form where the pocket is closed ( Neupert and Brunner 2002 



Neupert and Herrmann 2007). There are two differing opinions about the 



role of mtHsp70 in translocation. Either mtHsp70 pulls actively the substrate 
through the nanopore, or it only prevents backsliding of the substrate, leading 



to a passive mechanism (Click 1995 Neupert and Brunner 2002 Neupert 



and Herrmann[ |2007[). Such a ratcheting mechanism was first introduced by 



Schneider et al (1994) and Okamoto et al (2002). The best argument for active 



pulling is that the protein which is translocated must be unfolded outside the 



mitochondria, so translocation occurs against a force (Pfanner and Truscott 



2002). However, experimental evidence is still lacking and binding of mtHsp70 



to the substrate can explain observations ( Neupert and Herrmann 2007 ) 



Many proteins are translocated into the endoplasmatic reticulum, already 
during translation. Moreover, post-translational translocation has been sug- 
gested to occur through a passive ratcheting mechanism in eukaryotic cells by 



Matlack et al (1999). Here, the nanopore is formed by the Sec63 complex and 



the ratcheting molecules are called BiP which is a member of the Hsp70 family. 
BiP exists in an ATP-bound state with an open binding pocket. This is acti- 
vated by an interaction with the J-domain of Sec63, which can also close the 
binding pocket after BiP is bound to the substrate. Interestingly, BiP binds 
unspecifically to the substrate and therefore is able to mediate translocation 



of many different proteins (Rapoport 2007). For protein translocation into 
chloroplasts in plants, ratcheting mechanisms have not yet been considered 



yet (Strittmatter et al 2010 Li and Chiu 20101 



Quantitative descriptions of cellular ratcheting mechanisms began with pi- 
oneering work of Simon et al (1992) and Peskin et al ( |1993 ). In their model, 
ratcheting molecules can bind to the protein at the nanopore and hinder dif- 
fusion out of the nanopore completely (perfect ratchet) or only with a prob- 



ability p < 1 (imperfect ratchet). The assumption that ratcheting molecules 
can only bind directly at the pore is motivated by translocation into the endo- 
plasmatic reticulum where the ratcheting molecule, BiP, is known to interact 
with Sec61 which builds the pore. However, the interaction with the nanopore 
can only mean that concentration of ratcheting molecules is higher near the 



pore. (See |Griesemer et al (2007), for a mathematical model of the interaction 
of the ratcheting molecule with the nanopore.) We take the same approach 
as [D'Orsogna et al ( 2007 ) and assume that ratcheting molecules may bind 
uniformly along the protein. This assumes that the concentration of ratchet- 
ing molecules is approximately constant between the pore and the first bound 
molecules. 

Most quantitative descriptions of the Brownian ratchet assume a finite 



length of ratcheting molecules (e.g. Zandi et al 2003 D'Orsogna et al 20071 



Ambjornsson et al ( 1992 ) describe a car parking effect arising by the distribu- 



tion of ratcheting molecules along the substrate. They assume a fast binding 
and unbinding of ratcheting molecules, leading to an effective probability any 
position of the protein is bound. The case of fast binding is studied in |Bud- 



hiraja and Fricks ( i2006) , who derive a law of large numbers and a central 
limit theorem in the case that Brownian motion has a drift. Effects of bind- 
ing affinities which depend on the protein sequence have been put forward in 
[Abdolvahab et al| ( [20M| |20TT| . 



The goal of this paper is to introduce a realistic model for the ratcheting 
mechanism for protein translocation in mitochondria and the endoplasmatic 
reticulum as described above. .Elston ( |2002 ) tried to decide if translocation 
happens actively (by pulling the protein) or passively (by a ratcheting mecha- 
nism) using likelihood ratio tests. However, after fitting the model parameters 
using least squares, their statistical tests were not significant and both models 
fitted experimental data well. New molecular techniques have to give rise to 
better data to decide between active an passive transport. In the present pa- 
per, we obtain predictions for the speed of translocation which will shed light 
on empirical studies in the future. 

Our model of the Brownian ratchet is similar to the ratchet studied in 



Liebermeister et al (2001), but is continuous, assumes free diffusion of the 



protein and binding of ratcheting molecules anywhere along the protein. Since 
ratcheting molecules can dissociate from the protein, we call the process a 
broken Brownian ratchet. We give two models, termed Model I and Model II. 
For Model I, which takes into account every ratcheting molecule at all times. 



we can only show that the speed of translocation is positive (Theorem 2.1 1. 
Model II serves as an approximation. When a ratcheting molecule dissociates, 
we assume that binding and unbinding of ratcheting molecules is in equilib- 
rium. This model is easier to analyze but has about the same behavior for the 
speed of translocation (see the simulation results in Figure l3|. For Model II, 
we can show a law of large numbers giving the speed of transport in terms of 
a solution of a differential equation. Theorem |2.2[ 



2 The model 

The model we study extends the Brownian ratchet introduced in [Depper- 



schmidt and Pfaffelhuber (2010). Namely, protein translocation across a mem- 



brane is based on the following assumptions: 

(i) The protein moves in and out with the same probability. 

(ii) The protein movement is reflected at binding ratcheting molecules, 
(iii) The protein is infinitely long, 
(iv) The ratcheting molecules are infinitely small. 

(v) Ratcheting molecules may bind to the protein at a continuum of sites, 
(vi) The ratcheting molecules may dissociate from the protein. 



These assumptions are exactly the same as in Depperschmidt and Pfaffelhuber 



(2010) except for the last. They assumed that the dissociation rate of the 
ratcheting molecules from the protein is much smaller than their binding rate 
to the protein, leading to effectively no dissociation of ratcheting molecules. 

First, we translate the above set of assumptions into a mathematical model. 
By Xt, we denote the total length of the protein which already crossed the 
membrane between times and t. In addition, Rt < Xt is the largest distance 
of Xt to a ratcheting molecule. We assume that the protein has already moved 
across the membrane by time 0, such that ratcheting molecules may bind 
between Xt and — oo at constant rate. Consider the pair {Xt, Rt)t>o, where 
{Xt)t>o is a Brownian motion, started in Xq = x, reflected at Rt := swpTZt 
and TZt ^ (— oo, Xt] is the random set of all bound ratcheting molecules at 
time t. The dynamics of TZt is as follows: starting in TZq = {0}, a point at dx 
is added at rate j'i^xe{-oD,Xt]dx (i.e. a transition from TZt- to TZt '■— TZt- U {x} 
occurs). In addition, every x € TZt is deleted at rate S (i.e. a transition from 
TZt- to TZt '■— Ti-t-\{x} occurs at rate Stxi^Tit-)- Note that the last mechanism 
with rate S models the dissociation of ratcheting molecules from the protein 
(see (vi) above). We rely on the following graphical construction, which is 
illustrated in Figure [T] 

Definition 2.1 (7/^-broken Brownian ratchet, Model I) Let Af'^'^ 
be a Poisson point process on [0, 00) x M x [0, cx)) with intensity measure 
"/Se~^^X^{dT,dr,dz), conditioned on A/'^^*({0} x {0} x [0,oo)) — 1, where 
7 > 0, ^ > and X^ denotes the Lebesgue measure on [0,oo) x M x [0,oo). 
Moreover, J\f'' is the projection of M^- on the first two coordinates. Let 
{B^)t>o, {Bl)t>Q, ■ ■ ■ be a sequence of independent Brownian motions, in- 
dependent of JV' , with Bq = Bq = ■ ■ ■ — 0. We define recursively 
times to,ti,... when the reflection boundary changes as well as triples 
(to, ro, zq), (Ti,ri, 2:1), (r2,r2, Z2), • • • & N'^'^ such that r„ is the reflection 
boundary in the interval [i„,i„+i). Define to = tq = r^ = 0, zq :— z if 
(0, 0, z) £ Af^'^ and for xq > set Xt := |a;o + i3°| for to <t <ti, and the set 
of new reflection boundaries above r^ before time t„ + z^, 

r„ := {(r, r, z) G N''-^ ■ t,, < t < r„ + z„, r„ < r < Xr}, 
inf Fn '■— (t, r, z) if (t, r, z) G B^ and r = inf 7ri/^„, 



where tti is the projection on the first coordinate. Now, we set recursively for 
n = 0,l,2,... 

tn+i ■= (t„ + Zn) A inf 7rir„, 

'mir„, i/r„^0, 

{Tn+i,rn+i,Zn+i) := ^ {T,r,z) if Fn =0 and 

r = niax{r < Xt : t < Tn + Zn < t + z} . 

(2.2) 



Finally, 

Rt ■= Tn, 

Xt '■= Rt + \B^-t„ + Xt„- — Rt 



fortn<t<tn+i. (2.3) 



We refer to the process {Xt, Rt)t>o o-s the 7/5-brokcn Brownian ratchet, 
Model I with initial value {Xo,Ro) = (xo,0). 

Remark 2.2 (Interpretation) In the definition above, {T,r,z) S Af"*'^ rep- 
resents a ratcheting molecule which binds at time r at position r to the protein 
and stays bound up to time r + z. We start at time i = with one bound 
molecule at position 0. The variables ti,t2, ■ ■ ■ represent times when the re- 
flection boundary changes. Note that Rt, at any time t, is the position of the 
ratcheting molecule which is bound (i.e. t < t < t + z) and closest to Xt- 
By time t, let us denote by (r, r, z) € J^f-'^ the active point ii Rt — r (and 
recall that all points in M'^'^ have all coordinates diff'erent, almost surely). 
By definition, i„ is the time when the active point changes for the nth time. 
There are two possibilities when the active point {T,r,z) changes. First, a 
new Poisson point can fall between r and Xt (as described by the first line 



in (2.2)). Second, when t = t + z, the current active point dissociates, and the 



next active point is the one closest to Xt (given by the second line in (2.2)). 
Last, note that Xt^ = Xt^- by construction, for all n, meaning that {Xt)t>o 
is continuous and note that \Bf_^ +Xt^^ — Rt\ is a Brownian motion, started 
at Xt^- — Rt, reflected at 0, such that {Xt)t>o is reflected at Rt at all times. 
The graphical construction of {Xt)t>o is done step by step between jump 
times of the reflection boundary, {Rt)t>o, i-e. when the active point changes. 
This works since, after any finite time t there may only be a finite number of 
jumps of the reflection boundary so that the construction works for any large 
time. 

For Model I we show the following result. 

Theorem 2.1 (Speed of the broken Brownian ratchet, Model I) Let 

(Xt, Rt)t>o ^6 0, broken Brownian ratchet, Model I. Then, there are constants 
< c < C < oo, where c depends on 7 and S and C only depends on 7 such 
that 

c < a^ g := liminf — <a^s-= limsup — < C 

t-foo t ' t^oo t 

almost surely. Moreover, a^g — "f^'^aigi^2/3 anda^,s — j^'^aig/.y2/3 . 




^0 = To 



Fig. 1 One realization of the graphical construction of Model I. The path of {Xt)t>o is 
reflected at the boundary (-Rt)t>o- 



Remark 2.3 (Interpretation) The result says that the speed of the broken 
Brownian ratchet, Model I, (if it exists as a-y^s '■— Ql-j s — ^i,s) is positive, no 
matter how large 5 is. In addition, it scales with 7^' ^ like a^^s 



7^/^01,5/72/3. 



Although the speed of the broken Brownian ratchet is positive by Theorem |2.1[ 
we aim at a more complete picture. However, the difficulty in the analysis 
of Model I is that it is not local in the sense that a single Poisson point 
can be active more than once. In particular, dependencies between Poisson 
points arise as time evolves. For this reason, we consider a second model with 
similar properties as Model I. Briefly speaking, we introduce a Model II by 
assuming that the possible reflection boundaries below the currently active 
one are always in their equilibrium. 



Remark 2.4 (Motivation) Consider Model I and assume {T,r,z) is the ac- 
tive point by time t. Consider a time t' > t when the active point changes and 
i' = T + z. As an approximation, we assume that the set 

{r' : (r', r' , z') G Af'^ for some r', z' G [0, 00), t' < t' < r' + z'} 

is in its stationary distribution which arises for z — >■ 00, i.e. if the last active 
point stays for a long time. This stationary distribution is readily computed. 
It is clear that only points r' <r are possible. Moreover, the probability that 
a point in dr' is in the set equals 



76 



-Sz 



a/ 

dzdr' = —dr' . 




Fig. 2 One realization of the graphical construction of Model II. The path of {Xt)t>o i^ 
reflected at the boundary {Rt)t>0- Random variables Zq, Z\, ... are independent and Exp((5)- 
distributed, while Eq,E\, ... are independent and Exp(7/5)-distributed. 



So, the set is distributed according to a Poisson process with intensity measure 
^l(_oo,r](^')'^('^^')- ^^ ^ consequence, if the active point vanishes, the reflec- 
tion boundary jumps down an exponentiaUy distributed amount with rate 
7/(5 in the approximate modeL This leads to the following Model II, which is 
illustrated in Figure [2j 

Definition 2.5 (7/(5-Broken Brownian ratchet, Model II) Let M^ he a 

Poisson point process on [0, cx)) x M with intensity measure ^\^{dt.,dx), con- 
ditioned on H'~' {{0} X {0}) = 1, where 7 > and A^ denotes the Lebesgue 
measure on [0,oo) x M. Moreover, Zq,Zi,... are independent exponentially 
distributed with parameter 5 > and Eq,Ei,... are independent and ex- 
ponentially distributed with parameter ^/S (they are not needed if 6 = 0). 
Let [B'Dtx), {Bl)t>Oj ... be a sequence of independent Brownian motions with 
B° = Bl = • • • = 0. Define fo = r = and for Xq > set Xt := |a;o + B°\ for 
to ^ t < ti, o,i^d 



rn :- {(r, r)eW ■.tn<T< t„ + Z„, r^<r< X^}, 
inf r„ := {T,r) if {T,r) E r„ and t = miTTiFn, 



(2.4) 



where tti is the projection on the first coordinate. Now, we set recursively for 
n = 0,l,2,... 



n+l •- 



'n+1 




Zn) Ainf TTiTn, 

if rn^% and (tn+i-.r) 

ifTn^t 



inf Fn 



(2.5) 



Finally, 



Rt := Tn-, 



fort„<t<t,,+i. (2.6) 

t :— -tit -i- l^t-t ^ ^tr. 



Xt :=a + |R^- r +Xr 



We refer to the process {Xt, Rt)t>Q os the ^ / S-hrokcn Brownian ratchet, Model II 
with initial value {Xq,Ro) — (a;o,0). 

Theorem 2.2 (Speed of the broken Brownian ratchet, Model II) Let 

{Xt, Rt)t>o be a ^/S-broken Brownian ratchet. Model II. Then, 

t^oo t 2A{0) 

almost surely, where A is the first coordinate of a solution of the system 

A"{z) = -25B{z) + 27zA(z), 

B'{z) = -A'iz) 1b{z) ^^-^^ 

such that j4(0) = 1/2, i?(0) = and A is strictly decreasing with A{z) — >■ as 
z — > oo. 

Remark 2.6 (Simulations, uniqueness of A, and the case J = 0) 

1. Since Model II is only a convenient approximation of Model I, we use 
simulations to see differences in the speed of Model I and II. By scaling 



properties of both models (see Propositions 3.1 and 4.1 ), we require simu- 
lations only for a single parameter 7. As can be seen in Figure [3J the speed 
of both models is almost the same. The fact that Model II is faster for 
low 5 can be explained: we assume that the number of possible reflection 
boundaries below the currently active one is in equilibrium, which means 
that these are more than for Model I, where the equilibrium is not yet 
attained. Since more reflection boundaries mean that the broken Brownian 
ratchet moves faster, the speed of Model II is higher. For high values of 
5, Model I shows a higher speed in our simulations. The reason is that we 
fixed the reflection boundary at in our simulations of Model I, and hence 
the simulations overstimate the speed of protein translocation. 
For a proper use of Theorem |2.2| in Figure [3) we need to solve the sys- 



tem (2.7 1. We searched numerically for a solution by trying out various val- 
ues for j4'(0) and determined which value approximately leads to A{z) — )■ 
as z ^- 00. In our numerical analysis, we found only a single value of A' {Q) 
with this property, for all 5. Hence, we strongly conjecture there is a unique 



solution of the system ( |2?7| satisfying A(0) = 1/2, B{Q) = and A{z) -^ 
as z — )■ 00. 



3. In the case 6 = 0, the first equation in (2.71 reads A"{z) — 2^zA{z) = 



The only solution with the required boundary values is given by z h-> 
Ai{{2jy /^ z) / {2Ai{0)) , where Ai is the Airy function (i.e. solution of m"(z)~ 
zu{z) = 0) going to as z — >■ cx). By well-known properties of Ai (see e.g. 



Abramowitz and Stegun 1972 ) , the speed is thus given by 



lim — 

i— ^oo t 



Ai'jO) 
2Ai{0) 



r(2/3) /37\i/3 



(^^)''' = 7mu 



) « 0.36 • (27)l/^ (2. 



a result known from Theorem 1 in Depperschmidt and PfafFelhuber (20101. 



Theorem 2 
Model I 
Model 11 



Fig. 3 A simulation study for the speed of the broken Brownian ratchet, Model I and 
Model II. The parameter 7 = = is fixed. See also Remark 
represents the speed in the case 5 = from i2A 



The horizontal line at 0.36 



3 Proofs: Model I 



In this section, we prove Theorem |2.1[ Three ingredients are needed: First, a 
scaling property for the 7/(5-broken Brownian ratchet. Model I, is obtained in 
Section [XT] Second, an upper bound for the speed can trivially be obtained by 
setting (5 = and using results from Depperschmidt and PfafFelhuber (20101. 



Third, a lower bound is obtained by thinning Model I, which we do in Sec- 
tion |3.2| We conclude the proof in Section |3.3| 



3.1 Scaling property 



Here we obtain a scaling property of the j / S-hioken Brownian ratchet, Model I, 
which is based on a rescaling of space and time in Definition |2.1| 



10 



Proposition 3.1 (Scaling property, Model I) Let {Xf' , R]' )t>Q be the 
j/S-broken Brownian ratchet, Model I, with 7 > 0, (5 > and initial value 
(0,0). Then 

where — denotes equality in distribution. 

Proo f Using the notation from Definition |2.f[ we need to understand what 



(3.f) means for the Poisson process Af'^'^ and for the Brownian motions 
{B^)t>o, {Bl)t>Q, ■ ■ . We use the hnear rescahng of time and space 

^^' j(i,x,z) ^{jy%j^/^xn^/^z). 

and obtain g^iAf^^) = j\fhS/7^/\ Moreover, we have that ycB^)t>o = 
(5^4)00 foi' aU c > and k = 0, 1,2, . . . and so, if g^ : R"^ ^ K^ consists 
of the first two coordinates of g-y, 

g,{{t,B^),>o) = ({l'^\l'/'B^)\ 

"x / X (3-2) 

for k ~ 0, 1, 2, . . . where s :— 7^/'^i. Hence, the distribution of the Brownian 
motions are unaffected by the rescahng. We have shown that 

and the assertion follows. D 



3.2 Lower bound 

Consider the possible reflection boundaries [r, r+z) x {r} for all (r, r, z) S N^'^ 



in Definition 2.1 Clearly, the speed of the 7/(5-ratchet decreases if we take less 
possible refiection boundaries, leading to a lower bound for the speed of the 
broken ratchet. We will use the following thinned version of Model I. Note 



that the formal difference to Model I is the restriction supi7„_|_i < r in (3.4 1. 
An illustration can be found in Figure |4j 

Definition 3.2 (7/(5-broken Brownian ratchet, thinned Model I) 

Consider the same probability space as in Definition \2.1\ with the 
same J^'*'^ ,J^'' ,{Bf)t>o,{B})t>o,{B^)t>o, ■ ■ ■ Again, we define recursively 
times t(),ti,... when the reflection boundary changes as well as triples 
(t^O: H), ^0)1 (^1, ?^i, 2^1), (t^2, ?^2, ^2), ■ ■ ■ € Af^ such that r„ is the reflection 



11 

boundary in the interval [t„,t„_|_i). Define Uq :— 0, to ^ tq = rQ = 0, zq := z 
if (0, 0, z) e Af"''^ and for xq > set Xt := |a;o + 5°| for to <t <ti, and 

f„ := {(t, r, z) e TV"'''^ : i~„ < t < t„ + z„, f„ < r < X^} 
inf Fn := (r, r, z) if (r, r, z) € F„ and r = inf iriFn 

where tti is t/ie projection on the first coordinate. Now, we set recursively for 
n = 0,1,2,... 

tn+i ■■= (t„ + Zn) A inf 7ri_r„, 
^„+i := Sn U {inf{t e [f„, f„+i) : Xt = Rf}} 

iTn+l,rn+l,Zn+l) 

(t, r, z) if Fn — and 

r = max{r < X^ : sup Sn+i ^ t <Tn-\-Zn < t + z}. 

(3.4) 
fma/Zj/, 

^ ^ , ^ ^ Wori„ < i< t„+i. 3.5 

Xr.^Rt + \Bl^^+X^^^_^Rt\y "- 

We re/er to the process {Xt, Rt)t>o o,s the 7/(5-broken Brownian ratchet, 
thinned Model I with initial value {Xo,Ro) = (a;o,0). 

Remark 3.3 (Renewal times for the thinned Model I) Again, 
^0jiij^2, ■ • • are times when the reflection boundary changes for the thinned 
Model I. Moreover, the set Sn consists of all first times after changing the re- 
flection boundary, when Xt — Rt, up to i„. Clearly, the additional restriction 



supi7„+i < r in the definition of (i„+i,T„+i, z„+i) in (3.4) leads to taking 
less Poisson points and hence Rt < Rt and Xt < Xt for all t > 0, almost 
surely. Moreover, the times Z'oo := Un=o ^" ^^^ renewal times for the process 
{Xt — Rt)t>()- We have Xj — i?t = for t e X'oo and after time t, any reflection 
boundary r with (t, r, z) E JV^'^ must fulfill r > i in the thinned Model I. 
In other words, the possible Poisson points for the reflection boundary be- 
fore and after time t are distinct, so the distributions of {Xg — Rs)o<s<t and 
{Xs — Rs)s>t are independent. 

Using the renewal structure of the thinned Model I allows us to define a cu- 
mulative process. 

Lemma 3.4 (Renewal structure) Let {Xt, Rt)t>o be a broken Brownian 
ratchet, thinned Model I, and Sao — {cto,'?i, . . . }. Then, (fT„ — ct„_i, Xj^ — 
^CT„_i)ii=i,2,... ?s a sequence of bivariate iid random variables. 
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fTi ^2 h = T3 



Fig. 4 One realization of the graphical construction of thinned Model I. In contrast to 
Model I, there is the set of times Eao = {S'o i i?! i i?2 1 ■••} which serve as renewal points for 
the process. For example, note that the reflection boundary between 42 and (3 is lower than 
in Figure^ 



Proof By the construction of the thinned Model I, possible reflection 
boundaries r with (r, r, z) G A/''*'' are only used in one interval 
[ctqj'^i): [<?i,(T2)i ['?2j<?3), . . . Hencc, the Poisson points used for constructing 
((T„ — (T„-i, Xct„ — ^ff„_i)n>i arise independently, which shows the result. D 

Remark 3.5 (The thinned Model I as a cumulative process) 

Let (X(,i?t)t>o be the 7/(5-Brownian ratchet, thinned Model I. For Eao = 
{CTo,f?i,---}, set 



Mf, 



miuH n : a. 



>o, 



Sn '■— 2_^ \-^'Si 



^^^-l)l 



(3.6) 



so that we have Xt — SMt + ^t , where 



At := X: 



Co 



- Xt — Xg 



(3.7) 



According to the definition of cumulative processes in Roginsky (1994), the 



thinned broken Brownian ratchet A" is a type A cumulative process with re- 



Xg 



are 



mainder At- Note that finite first moments of cti — ctq and X-g^ 
sufficient for the strong law of large numbers for Smi ■ We will see that a law of 
large numbers holds also for (Xt)t>Oj as we will show in Proposition 
At/t converges to almost surely as t — ^ 00. 



3.8 



that 
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Lemma 3.6 (Finite first moments of cr) Consider a thinned Model I, 
started with Xq — x > and 

do ■■= inf{i >0:Xt= Rt}, (3.8) 

that is, (To *s the first time when then broken Brownian ratchet hits the moving 
reflection boundary. Then, using Ex['] '-^ ]E[-|Xo = x], 

supExiao] < oo. 

X 

Proof Set Ht :— Xt^Rt and note that H — {Ht)t>o locally behaves like Brow- 
nian motion, and has jump discontinuities which, by time i, either increase or 
decrease Ht-. First, Ht- jumps to h < Ht- at rate ''ydh by occurrence of a new 
Poisson point in the graphical construction at rate 7. Second, it jumps to some 
h > Ht- at rate S since t = t + z for the active Poisson point (r, r, z) by time 
t—. Note that ctq := inf{t > : Ht — 0}. We couple the process {Ht)t>o to a 
process {It)t>o with Ht < It for aU < i < do, define ct/ := inf{i > : It — 0} 
and show that sup^, E^; [cr/] < cx). 

The process I = {It)t>o has the following dynamics: it behaves locally like 
the same Brownian motion as "H, but starts in /p = 00. It jumps to 1 at time t 
in the following cases: (i) if H jumps from Ht- > 1 to Ht < I (ii) if H jumps 
from Ht- < 1 to Ht < Ht- (iii) at an independent rate 7(1 — Ht)^. In total 
this gives a jump rate of 7 to jump to 1 and Ht < /* at such a jump time. In 
addition, I jumps to 00 if H jumps from Ht- to Ht > Ht-, which occurs at 
rate S. 

Consider an exponentially distributed random variable T with rate "f + 5 
and a Brownian motion, starting in 1 at time 0. Let as be its hitting time of 
and set p := V{aB <T)>0. Note that c/ is given as follows: A Poisson point 
at rate 7 has to occur at time s with the property that the Brownian motion, 
starting in 1 at time s hits before any of the other events, bringing I back 
to 1 or infinity, occurs. We obtain 

E.,[ai] = — + EiasWB < T] = — + ^ ^' ^ < - ( - + ——, < 00, 

P7 P7 p p\j J + 0/ 

independently of x. The result follows since do < <7i, almost surely. D 

Lemma 3.7 (Increments of {Xt)t>o) For the j/S-broken Brownian 
ratchet, thinned Model I, started in x > 0, we have 

(t) E,[X|J < 00, ^ 

(ii) <Ex[Xs^ -Xffo] < 00 7 
(iii) < E^ [ai — ao] < 00 . 

Proof For (i), note that I^Sol — l^pol- Using the second Wald identity and 



Lemma 13.61 we obtain 

E,[X|j=E,[Bij=E,[ao]<oo, 
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which shows (i). For (ii) and (iii), we can assume that a; = and ao — 
without loss of generahty. Consider the time t of the first jmnp of {Rt)t>o- Of 
course, we have X^- > almost surely and therefore < E[^f] < oo. Moreover, 



we have that (^s^At ~ ^t 



is a uniformly integrable martingale. Indeed, 



again using the second Wald identity and Lemma |3.6[ 

supE^iXs^M - Xi~)^] =supE4aiAt-t] =E^[di-t] < oo. 
t>T t>t 

By the Optional Stopping Theorem we obtain that 

EoiXsJ = Eo[Xy] + Eo[Xs, - X^] = Eo[Xj.] 

and (ii) follows since < E[Xijr] < oo. For (iii), rewrite cti — (di ~ t) +t. We 
know that t is bounded from above by the killing time r of a Brownian motion. 



if it is killed at rate j\Bt\ at time t. In Depperschmidt and Pfaffelhuber (2010 



Lemma 5.4) it was shown that Eq[t] < oo. Hence, by Lemma 3.6 

Eo[di]=Eo[di-t\+Eo[t] 



Eo 



E 



Xt-R, 



[f^oj 



-Eo[i] <supE^[CTo]+Eo[r] < oo, 

a:>0 



which finishes the proof. 



Proposition 3.8 Recall from Remark 3.5 the j/6-broken Brownian ratchet 



as a cumulative process with remainder {At)t>o from (|3.6[). Then 



lim — 

i— >-C30 t 



almost surely. 



Proof We proceed as in Lemma 8 in Smith ( 1955 1. Define for n = 1, 2, . . . 

y„= sup \Xt-Xsj 

te[CT„-i,CT„] 
and note that Yi, I2: ■ • • are independent and identically distributed. We write. 



using Lemma 3.7 i) and the law of large numbers for {Mt)t>o 

\Xt - Xs 



hm ^ < hm ^ 

i— >-oo t i— >oo t 



lini 



t 



"*' < lim ^. lim ^ 

t— 4-00 t i— >-oo A^i 



1 ,. r„ 

= — 7^ ;^T ■ hm — . 

E[o'i — (ToJ n^oo n 

Since Yi < X-g^ + sup^Q Xt^^i ■, by Doob's inequality, 

E[Y^] < 2E[X|J + 2E[supXf^-J < 2E[X|J + 8E[X|J < 00, 



(3.9) 



it follows that Y„/n — > almost surely, as n — >■ 00 by Lemma 7 in Smith 
(|1955|). Now, the result is implied by (pjl). D 
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3.3 Proof of Theorem [2H 

The proof consists of three steps. We derive the lower and upper bound for 
the speed of the broken ratchet. Then, we prove the scahng property. 

Step 1: Lower bound for the speed: The thinned Model I of the 7/(5- 

broken Brownian ratchet uses less possible reflection boundaries, leading to a 
smaller speed. Hence, if we can show that it still has positive speed we obtain 



that Model I has positive speed as well, as claimed in Theorem 2.1 

In order to prove that thinned Model I has a positive speed, we use the 
law of large numbers for the cumulative process as introduced in Remark |3.5| 
By the law of large numbers for (5'„)„=i_2,... and (Mt)t>o, Lemma 3.7 and 
Proposition |3.8[ we have 



Xt 5'a/j Mt At t-foo ^[Xg — Xs ] 

— = -TT > ^r^ ^ , =: c almost surely. 

t Alt t t E[cti - CTo] 

Step 2: Upper bound for the speed: The speed of the -f/S-hroken Brow- 
nian ratchet is decreasing in S. In particular, taking 5 = gives an upper 
bound for the speed. Moreover this bound is independent of S. We find from 



Depperschmidt and Pfaffelhuber (2010) that 



— - — < — > C almost surely. 

t - t ^ 

for some < C < oo. 

Step 3: Scaling property: As a consequence of Proposition |3 . 1 1 we have 



liminf— !^ =liminf7 3 = liminf 73 , 

t— Xxj t t— foo t t— >-co ^^l-^t 



.l;5/7^/^ 



= lim inf 7 3 

t-»-oo t 

where lim inf can as well be replaced by lim sup. 



4 Proofs: Model II 



We start in Section [4~T| Proposition 4.1 with showing the same scaling prop- 



erty as given in Proposition |3. l| for Model I. The rest of the section is concerned 
with concrete calculations. Model II is linked to killed Brownian motion in Sec- 
tion 4.2 We interpret Model II as a (random) additive functional of a Markov 



chain, where the related Markov chain is defined in Definition HTS^ and is shown 



to have a unique equilibrium in Propositions 4.7 and |4.8| in Section |4.3| For 
this equilibrium, it is possible to compute increments analytically, which is 
carried out in the proof of Theorem |2.2| in Section |4.4[ 
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4.1 Scaling property 



The following result is the analogous result of Proposition 3.1 for Model II 



Proposition 4.1 (Scaling property, Model II) Let {X^' , R]' )t>o be the 
^/S-broken Brownian ratchet, Model II, with 7 > 0,(5 > and initial value 
(0,0). Then, the following holds: 

{^t l-Kj jt>0 - 7 ^ (^A^2/3t '^72/3f ) ■ (4.1) 



Proof We use supercripts 7, 5 in order to emphasize dependency on the pa- 
rameters and start with the {X]' , i?7' )t>o-ratchet. Recall the notation from 

the proof of Proposition 



3.1 



as well as gjN "') = N^. Note that the laie-S ran- 



dom variables Z^^^Z^-,- ■ ■ in Definition 2.5 only affect the i-direction, while the 
random variables El' , £"2 ' i -^3 ' j • ■ ■ ^-t rate 7/(5 affect the x-direction. Since, 



7& T^^-M ± i^ih^'^ pi;5/7^/'- 



gj{Zl[,E^^' ) = (Zn tEu ), and since the same rescaling as in (3.2) 

holds in Model II, we have shown that 






[^t ^^t )t>a — \^^2/'n '^72/34 Jt>o 
and the assertion follows. 



4.2 Connections to killed Brownian motion 

The Brownian motion which drives the broken Brownian ratchet. Model II, is 
reflected at the same boundary until Rt changes. We interpret the time when 
this happens as a killing time of the Brownian motion. Due to the scaling 
property of the Brownian ratchet it would be possible to carry out the com- 
putations in the notationally convenient case 7=1 and deduce the general 
result from this particular one. However, we feel that an approach explicitly 
allowing for all values of 7 is more transparent. 

Definition 4.2 (Killed reflected Brownian motion) Let {Bt)t>Q denote 
Brownian motion started in x > and consider the stopping time 77 defined by 



inf{t>0: /" {-f\Bs\+S)ds>^y 



7? := inf |t > : / (7|S,| + 5) ds > ^j, (4.2) 

where ^ is an independent exponentially distributed random variable of rate 1 . 
Define B := (i?t)t>o by Bt := \Bt\ for^O <t<V and Bt = A for t > r], 
where A ^R is the cemetery state. Then B^ = B is reflected Brownian motion 
killed at rate ^Bt+5 at time t. Denote the probability measure of the Brownian 
motion started in x by V^ and write E^, for the respective expectation. 
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The second order operator A associated to the killed Brownian motion defined 
above acts on C^ functions / : [0, oo) — >■ M satisfying /'(0+) = according to 



Afix) = -nx)-{jx + S)fix). 



(4.3) 



The diffusion process B is transient since it can be killed in each interval with 
positive probability. Let p(-;-,-) denote the corresponding transition density 



m 

by 



with respect to the speed measure which is given by m{dx) — 2dx (see Borodin 
and Salminen (2002[ p. 17)). In the same way as in the proof of Lemma 5.4. 



Depperschmidt and Pfaffelhuber (2010), the Green function of B^ is given 



G{x,y) 



p{t\x,y)dt, 



which is finite for each x, y > 0. It can be written in terms of solutions of the 
differential equation Au = which reads 



u"{x) - (27a; + 25)u{x) = 0. 
The following lemma is the basis for all calculations to follow. 



(4.4) 



Lemma 4.3 (Green function of killed Brownian motion) The Green 
function of the killed Brownian motion from Definition \4-S\ is given by 



(4.5) 



1 j -(/'(a;)0(y) forO<x <y, 
G(x, y) := — < 

w \^il){y)(l){x) forO<y<x, 



where 



w 



(27) 



1/3 



(4.6) 



with 



(t>{x) =Ai({2-fy/^x + 
ipix) = Bi({2"fy/^x + 



2^'H 



C = C{-i,5) 



Bi^(^ 
^*'( ^2/3 ) 



are two solutions of (4.4) 



CA 



^((27)' 



/3. 



(4.7) 
(4.8) 

(4.9) 



Proof In order to compute the Green function of the killed Brownian motion 



we need two particular solutions, say / and g, of the equation (4.4) on [0, 00) 
such that (see e.g. Borodin and Salminen 2002 p. 18-19) 
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(i) / is positive, strictly decreasing and f{x) — )• as a; — )• oo, 
(ii) g is positive, strictly increasing, 
(iii) g'{0) = (this is the condition for reflecting boundary). 



The functions x t-^ Ai 



((27)^ 



/'x+?^ 



and X i-^ Bi 



((27)^ 



/3.T + 21^ 



are two linearly independent solutions of (4.4 1. It is easy to check that the 
requirements (i)-(iii) are satisfied by functions f = (p and g — ip defined in 
(lijl) respectively (li^Sl. 



The Wronskian of ^ and (j) is given by 
w = '(/''(O)0(O) - V'(0)(/)'(0) = -V'(O)0'(O) 

= {2Jy/'{B^\'^)A^C^) - B^C^)A^'C^)) 



(27)1/3 (4.10) 



hence all results follow from (see e.g. Borodin and Salminen 2002 p. 18-19). 

D 

The density of the killing position with respect to Lebesgue measure is k{y) — 
2"fy + 26 and moreover, the killing position of the Brownian motion started in 
X has density 



f{y):=G{x,y)k{y) 



(4.11) 



with respect to Lebesgue measure (see |Borodin and Salminen 2002 p. 17 resp. 

p. 14). 

The next lemma will be used in the proof of Proposition |4.7| 



Lemma 4.4 (Bounds on increments of killed Brownian motion) For 

any x > 



'tm^<EAB.,^]<x+'tmm, 



w 



w 



(4.12) 



Proof Using the fact that ip and cj) solve (4.4) and then integration by parts 
we obtain 



EABr,-]= / yG{x,y)k{y)dy 



^^""^ f^y^{y){2^y + 2S)dy+^ I yij{y){2jy + 26) dy 



= t^ r ycl>"{y) dy+^ r yriy) dy 

W J^ W Jo 

= ^(-V^(x)0'(y) + 0(x)V.'(.)) + '^^'^M^ ^ X + -^(^MO) . 

W WW 

Now using the fact that (j) is strictly decreasing we obtain the upper bound 
in (4.12). Furthermore, due to ^'(0) = by the choice of C we have -0(0) = 
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-w/4i'{Q) = -w/v4i'(^^). As the function 



a; I— >■ X - 






is strictly increasing on [0,oo) its minimum is attained in a; = 0. This gives 
the lower bound in (4.12[). D 



4.3 The invariant distribution for increments at jump times 

One advantage of Model II as compared to Model I is that we can define 
a Markov chain which models time- and space-increments at jump times of 
the reflection boundary. We show that the Markov chain has an equilibrium 



(Proposition 4.7) which is unique (Proposition 4.8). 



Definition 4.5 (Markov chain at jump ti mes ) Consider the ^/6-hrok 



en 



Brownian ratchet {Xt,Rt)t>o from Definition 2.5 with Xq = x and Rq = 0. 
Define {Yn,Wn, r]n, )n>i by 

^n := ^f„ - fn, Wn := r„ - f„_i and ?/„ := t„ - i„_i. (4.13) 

Note that for any k, {Yn,Wn,r]n)n=k+i,k+2,... depends on (F„, W„, ry„)„^i^2,....fe 
only through Yfe. That is, {Yn,Wmrin,)n>i is a Markov chain. For any n>\ 
we denote the law of {Yn,Wn,rin)n>i by (P"). 

Remark 4.6 (Distribution of Yn and Wn) Let (£'„)„=o,i,... be iid expo- 
nentially distributed with parameter j/S and let (C/„)„=o,i,... be iid uniformly 



distributed on [0, 1]. Then, recalling Definition 2.5 



V — 

-'n+l — 



(^t„+i - ^n)(l - Un), if r„+i > f„, 
{X^^_^^ -rn + En), if r„+i < f„, 



and 



I -En, if r„+i < f„ 



d 



where = means equality in distribution. In particular, on the event {Pn+i > ?„} 
the random variables Wn+i and Yn+i have the same distribution. 

In the following two propositions we obtain existence and uniqueness of the 
invariant distribution of the Markov chain (Yn, Wn, ?7n)n>i- 

Proposition 4.7 (Existence of invariant distribution) For any x > the 

sequence of Cesdro averages (n^-^ X]fc=i ^x) >i converges along a subsequence 
weakly to an invariant distribution P^ of the Markov chain {Yn,Wn,iln)n>i- 



20 



Proof In the case 5 = the assertion is shown in Proposition 5.6 in |Dep- 



perschmidt and Pfaffelhuber (20101. Thus, we assume i5 > in the rest of 



the proof. We only need to show that the first moments of F„, W„ and rjn 
are bounded uniformly in n. This imphes that the sequence (-P")n>i is tight. 
Then the sequence of Cesaro averages [n~^ Sfc=i ^x)n>i ^^ ^^^'^ tight. That 
is, any subsequence contains another subsequence along which the Cesaro av- 
erages converge weakly. These weak hmits are invariant for the Markov chain 
(y„, Wn, rin)n>i- A proof of this fact in the continuous time case, that can be 
easily adapted to the discrete time case, can be found in Liggett ( 1985[ p. 11). 
The moments of each 77„ are bounded from above, since by construction, see 
(2.5), each r]n is bounded by an exponential random variable with parameter 



6. Furthermore by construction we have (see Remark |4.6[ ) 

Yn+l = (-'^tn+i ~ ^n)(l{f„ + i>f„}(l - Un) + l{f„ + i<f„}) + l{f „ + i<f „}-E'n 
< {Xi^^^ - f„)(l - t{f„^,>f„}Un) + En. 

We set Qn = a{{Xt,Rt)o<:t<i„) and c := (/)(0)V'(0)/w. By the strong Markov 
property we obtain 



where Bt is killed reflecting Brownian motion from Definition 4.2 and ry is its 
killing time. Furthermore 



E[{Xi^^^ - fn)t{f„^,>,^}\Gn] = EyJiXi^ - fo)t{,-,>ro}] = ^Y^ 



K 



B„ 



5h 



The function x M- x^ / (x + 5/7) is convex. Thus, by Jensen's inequality 



E[(X 



t-n + l 



^n)l{f,.+i>f„}|5n] > Ey„[-B,,-] 



EyjBr,- 



Ey„ [B„_] + (5/7 



Now using that the function x M- x/{x + (5/7) is increasing and (4.12) we 
obtain 



s/Y 



Altogether, setting 



q:=l- 



1 



e(0,i) 



2 c + (5/7 
we obtain 



Un 



En 



= E, 


E[(^t„+, - 


f„)(l-l{,„^,>,,.}f/„)lg„] 


<gE, 


Ey„[S^-] 


+ 6h 


<q{E, 


AYn]+C)^ 


S/l, 



(5/7 
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where we used (4.12 ) for the last inequaUty. Uniform boundedness of the first 
moments of Yn follows since the recursion x„+i = qx„ + d has a unique fixed 
point for |g| < 1 and any d. For Wn+i we have 

Wn+l = l{f„ + i>f„}C/n(-'^t„^i - fn) ~ l{r^^ j<r^}il^„. 

Now using that on {r„+i > f„}, Wn+i — |W„+i| and K„+i have the same 
distribution we obtain 

E.[|w^„+i|] = E,[i{,„^^>,„}y„+i + 1{^,„^^<,„}S„] < E,[y„+i] + s/j. 

Uniform boundedness of first moments of Wn follows from uniform bounded- 
ness of first moments of Yn. D 

Proposition 4.8 (Uniqueness of invariant distribution) Consider the 



Markov chain (Yn, Wn, '7n)n=i,2,... from Definition 4-5 on the basis of a broken 
Brownian ratchet, Model II started in x > 0. If an invariant distribution of 
the Markov chain exists, then it is unique. 

Proof In the proof we use a second graphical construction for the 7/(5-broken 
Brownian ratchet, Model II, {Xt, Rt)t>o, which is based on a single Brownian 
motion (-Bt)t>o with Bq = 0; see also Figure IslJA). 

We define the sequence of active points {Sn)n=o,i,..., the sequence of re- 
flection boundaries {Rn)n=o,i,... and a sequence of jump times {Kn)n=o,i,... 
as follows: Let (f„)„=o,i.... be iid exponentially distributed with parame- 
ter 1, {En)n=o.i,... iid exponentially distributed with parameter ^/S, and let 
(Un)n=o.i.... be iid uniformly distributed on (0, 1). We define kq — 0, Sq = x, 
i?o = and for n > 

Kn+i = inf |t > K„ : / {'j\Bs - 5„| + S) ds > S,n > 

is the time of the next event after k„ with rate 7|i?s — ^nl + (5 at time s. At 
time Kn+i, there are two possibilities: 

• If {^„ > S,„^ J, set 

{{Sn — Un{Sn — -Bk„+i), -Rn + Un{Sn — Bk,„+i)) 
with probabihty ^^^^^^"+lg^J'^^ , 
(Sn + En, Rn " En) with probability ^[^^^ ^^-s„\+s - 

• li{S„ <B^^^,}, set 

(Sn + C/„(-Bk„^i — Sn),Rn + UniB^^^-^ — Sn)) 

{Sn+i,Rn+i) = { with probability ""+' 



{Sn — En, Rn — En) with probability 



s 



22 



Thus, the active point jumps uniformly between the currently active point and 
the Brownian motion at rate 7 times their distance or it jumps an exponential 
distance away from the Brownian motion. In particular, 



Rn+i — 



Rn + CAi|-Bk„+i — Sn\ with probability 



7|S»„_^i-S„| 



^1 Un\ WiUll piuucxuiiitj ^|^^^^^_5,^|^^-, 



Rn — En with probability 



s 



i.e. the reflection boundaries either jump up at rate 7 times the distance of 
the Brownian motion to the active point or down at rate 5. 
Now set 



Rt = Rn, St = Sn, for t e [Hn, Kn+l) 



and 



Xt ^ Rt + \Bt - St\ fori>0. 
Then, {Xt)t>o is reflected at {Rt)t>o and is continuous, since for all n > 1, 

^K„+i— ^ Rn I |-^K„+i— Sn\, 

Rn + C^kI-Bk^^i- — Sn\ + |-B„,^_|_j — Sn ~ Un{Bi^^^-^ — Sn)\ 



^K„ + l 



with probability ''"^^ 



Rn — En + |-Bk„+i — ^nl + i?„ 

with probability -rg „ , , , 

In other words, the process {Xt, Rt)t>o is a j/S-hioken Brownian ratchet. 
Model II starting in (|x|,0). 

We now prove Proposition |4.8| using a coupling argument, illustrated in 
Figure [sIB). Recall that a coupling of two 7/(5-broken Brownian ratchets. 
Model II, is a process {X^^\R'i^\xP ,r['^^) such that (Xj^^^flf^) is a j/S- 
broken Brownian ratchet, i = 1, 2. It is straightforward, using the same Brow- 
nian motion, to construct such a coupling along the lines of the above con- 
struction, where {X^^\r[^^) starts in {x^^\0) and (xf\^P^) in {x^^\0), 
respectively, for x^^\x^^^ > 0. Then, we must show that both broken ratchets 
use the same active points after an almost surely flnite time. This sufSces, 
since from this time on the time- and space-increments at jump times are the 
same for both processes, and hence, both must converge to the same limit law 



(see e.g. [Lindvall] [1992| Thm. 21.12). 

In order to show that both broken ratchets use the same active points after 
finite time, note the following: if the active points of both ratchets are both 
above or both below the Brownian motion they have a chance to jump to the 
same point (see Figure IsJB)). Assume for definiteness that at time t we have 



Sf' < SI ' < Bt (as at time t in Figure^B)). Then at rate -i{Sl' -Si') only 
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the lower active point jumps towards the Bt and if it does, it jumps uniformly 
between S^ and S^ . Furthermore at rate ^{Bt — Sl ) both jump to the same 
uniform position between SI and Bt- From then on the active points evolve 
together and we refer to this time as the coupling time {Tc in Figure p[B)). 
After the coupling time both ratchets have the same jump times and the same 
increments. Since the Brownian motion will spend infinite amount of time 
below or above both active points the coupling time is almost surely finite. 

D 
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Fig. 5 (A) Alternative graphical construction of Model II from the proof of Proposition 
(B) Graphical construction of two coupled 7/5-broken Brownian ratchet, Model II, using 
the same Brownian motion, similar to (A). After time Tcoupi the active points and the jump 
times of active points are the same for both processes. 



4.4 Proof of Theorem 



The proof of Theorem [2^2] comes in three steps. First, we will connect the speed 
of the broken Brownian ratchet, Model II, to the Markov chain (F„, Wn, ??n)n>i 



from Definition 4.5 and show that the speed equals the ratio of expected 



space- and time-increments between jump times. In Step 2, we use the unique 
equilibrium of this Markov chain to obtain a fixed point equation for Y, which 
allows to study its density. In Step 3, we use insights into the density of Y in 
order to study expected increments between jump times. 



Step 1: Connection to the Markov chain (1^, Wn,rin): In the last sec- 
tion we have established the unique equilibrium distribution n for the Markov 
chain {Yn,Wn,'rin)n>i- Let {Y,W,r]) be distributed according to tt and recall 
that Y is the distance of the Brownian ratchet to its reflection boundary at 
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jump times, W is the amount by which the reflection boundary jumps and 77 is 
the time between jumps. Note that 77 equals the killing time of the Brownian 
motion as introduced in Definition |4.2[ if the Brownian motion is started in 
an independent copy of Y . 

The proof works along similar lines as in the proof of Theorem 1 in |Depper^ 
Schmidt and Pfaffelhuber (12010). Formally, Model II is a cumulative process 
similar to the thinned Model I; see Remark |3.5[ Renewal points are times 
when the Brownian ratchet hits the reflection boundary, given a jump of the 
reflection boundary occurred between renewal points. The law of large num- 
bers for cumulative processes establishes that Xt/t converges as i — ?> 00 to the 
ratio of expected space- and time increments between renewal times, almost 



surely; compare with (4.14) below. Then, proceeding exactly as in the proof of 

923) 



Theorem 1 in Depperschmidt and Pfaffelhuber (2010 



from the ratio limit theorem (see e.g. Revuz 1984 



P- 



it can be seen 
Thm. 6.6 on p. 231) that 



,. Xt E[W] 
hm — = -— r- r. 
t-)-oo t E[rj] 



(4.14) 



Step 2: A distributional identity and some consequences: The triple 
(y, W, rj) is distributed according to the equilibrium tt. Starting the Brownian 
ratchet in position Y, the reflection boundary at the next jump time either 
jumps up or down. Let U be uniformly distributed on [0, 1] and denote the 
length of the possible downward jump by Z, distributed as Exp ((5/7). By 
construction 



\^ — 11 R^ Tl 

-* -"^{junip up at r]}^r]— 



' -"^{jump down at 77} 



Bl 



(4.15) 



Conditioned on the event that the reflection boundary jumps at time 77 the 
probability of an upward jump is 



iB, 



Bl 



7<_+5 B^_ 



5h 



and for a downward jump 



5/1 



jB^_+S B};_+dh 
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and (4.11 1 imply 



Let us write fx for the density of a random variable X. Equations (4.15) 

friz) ^ I fvix) I /5. {u) (^^fu (-) - + -^^ 
fvix) / /bx iu) .r^dudx 



fz {z — u) \ du dx 



^^- u + (5/7 

+ r Iy{x) r f^A^)^l^^e-(^--^^/'dudx 
Jo Jo " u + 5/7 d 



/>00 /"OO 

= 2j fvix) / G{x,u) dudx 

'• ^-v ' 

= :A(z) 

pOO pZ 

+ 27/ fvix) G{x,u)e-^'-"^''/^dudx. 
Jo Jo 

> - -/ 

= :Biz) 

In particular we have 

/■OO /"OO 

/y(0) = 27^(0) = 27 / /y(x) / G{x,u)dudx. (4.16) 

Jo Jo 

We aim at showing that the pair (A, B) satisfies the differential equation ( |2.7[ ). 
For that we need to differentiate A and B. The derivatives of A and B are 
readily computed as 

1 



A'{z) 



and B'{z) = -A'{z) - ^B{z). 
5 



■0(2) / /y(a;)0(a;)da; + 0(z) / fy {x)^l^{x)dx 



7 



Thus, 



For A we have 
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/;.(^)^27(A'(z) + S'(z)) = -^B(z). 



^(^) = / fvix) I G{x,u)dudx, 
A'iz) = 
A"{z) =. - 



A"\z) 



1 

w 
1 

li; 

1 r 

w 



'4>{z) I fY{x)(l){x)dx + (t){z) fY{x)ifj{x)dx 



^'{z) fY{x)4>ix)dx + <j>'{z) fY{x)i,{x)dx 



(4.17) 
(4.18) 

(4.19) 



V/'(z) / fY{x)<t>{x)dx + 4>"{z) fY{x)^{x)dx 

'- J z Jo 

^fY{z){^'{z)cj,{z)-i^{z)<l,'{z)) 
/y(z) + (27Z + 2J)A'(z), 



26 



where we have used that 4> and tp are solutions of (4.4), the equation for A' 
and that the Wronskian does not depend on z. Writing now fy = 27(A + B) 
and using (4.17) we obtain 

A"'{z) = 2-f{A{z) + B(z)) + (27Z + 25)A'{z) 

= -2S{A'{z) + B'{z)) + 2^A{z) + {2^z + 25)A'{z) 
= -25B'{z) + {2-^zA{z))' . 



Integrating both sides we have 

A"{z) = ~25B{z) + 27zA(z). 



(4.20) 



The integration constant vanishes because -B(O) = and A"{Q) = 0, as can be 
seen from ( |4.19[ ) (recall that V'(0) = 0). Thus, [A, B) is a solution of ^7j\ . 

Step 3: Computation of ]E[VF] and E[?7]: In the case that the reflection 
boundary jumps up W and Y have the same distribution. If the reflection 
boundary jumps down W is distributed as —Z where Z is exponentially dis- 
tributed with parameter 7/(5. Thus, the density of the invariant distribution 
of W is given by 

Iw{z)^ I Iy{x) I f§. (u)( t{z>o}fu (-) - 

Jo Jo "- Vu + (5/7 ^-' \uJ u 

u + (5/7 



/2(-2)l{2<o} ) dudx 



= 27A(z)l{,>o}+y fvix) j f^^^_{u)-^^'^-e'^"t{,^o}dudx 

= 27A(z)l{,>o}+e-^/^l{,<o} / /y(x) / G(x,u)dudx 

Jo Jo 

= 27A(z)l{,>o} + 27e^''/^l{,<o}A(0). 

In order to compute E[W] we have to integrate the last equations. We use 
(|4.20|) and (|4.17|) and obtain 



E[W] = / 2jzA{z) - 2jze^^^/^A{0)dz 
Jo 



/>oo 

/ A"(z) + 26B{z)dz 
Jo 



2^ 

7 



A(0) 



(4.21) 






Moreover, by (4.16) we have 

EM = 2 



/y(x) / G(x,y)d2/da: = 2A(0). 
Jo 



(4.22) 
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Here the factor 2 shows up because we have to integrate the Green function 
with respect to the speed measure, which is m{dy) = 2dy in the case of 
Brownian motion. 



Now, Theorem 2.2 foUows from (4.14), (4.21) and (4.22), since multiplying 



A by a constant factor does not change the ratio —A'{0)/{2A{0)). 
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